Peculiar Velocity Decomposition, Redshift Space Distortion and Velocity 
Reconstruction in Redshift Surveys — I. The Methodology 
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Massive spectroscopic surveys will measure the redshift space distortion (RSD) induced by galaxy 
peculiar velocity to unprecedented accuracy and open a new era of precision RSD cosmology. We 
develop a new method to improve the RSD modeling and to carry out robust reconstruction of the 
3D large scale peculiar velocity through galaxy redshift surveys, in light of RSD. (1) We propose a 
mathematically unique and physically motivated decomposition of peculiar velocity into three eigen- 
components: an irrotational component completely correlated with the underlying density field (v^), 
an irrotational component uncorrelated with the density field iys) and a rotational (curl) component 
(v_b). The three components have different origins, different scale dependences and different impacts 
on RSD. (2) This decomposition has the potential to simplify and improve the RSD modeling. (I) 
vs damps the redshift space clustering. (II) vs causes an overall damping. But it also enhances 
the redshift space power spectrum P''{k,u) with a leading order contribution of directional 
dependence. Here, u = k^/k, k is the amplitude of the wavevector and kz is the component along the 
line of sight. (Ill) -vg is of the greatest importance in the RSD cosmology. However, even in the limit 
of V5 — >■ and vs — ^ 0, the leading order term of P''{k,u) is cx (I + fW{k)u'^)'^. It differs from the 
usual Kaiser formula by a window function W{k). Nonlinear evolution generically drives W{k) < 1. 
We hence identify a significant systematical error causing underestimation of the structure growth 
parameter / by as much as 0(10%) even at relatively large scale k = O.f/i/Mpc. (IV) The velocity 
decomposition reveals the three origins of the finger of God (FOG) effect and suggests to simplify 
and improve the modeling of FOG by treating the three components separately. (V) We derive a 
new formula for the redshift space power spectrum. Under the velocity decomposition scheme, all 
high order Gaussian corrections and non-Gaussian correction of order S"^ can be taken into account 
without introducing extra model uncertainties. (3) The velocity decomposition clarifies issues in 
peculiar velocity reconstruction through 3D galaxy distribution. We discuss two possible ways to 
fulfill the 3D va reconstruction. Both use the otherwise troublesome RSD in velocity reconstruction 
as a valuable source of information. Both have the advantage to render the reconstruction of a 
stochastic 3D field into the reconstruction of a deterministic window function W''{k,u) of limited 
degrees of freedom. Both can automatically and significantly alleviate the galaxy bias problem and, 
in the limit of a deterministic galaxy bias, completely overcome it. Paper I of this series of works 
lays out the methodology. Companion papers [1| will extensively evaluate its performance against 
N-body simulations. 

PACS numbers: 98.80.-k; 98.80.Es; 98. 80. Bp; 95.36.+X 



I. INTRODUCTION 



The observed galaxy clustering pattern in redshift 
space is modified by peculiar velocity of galaxies and 
shows characteristic anisotropies 043] • This redshift 
space distortion (RSD) effect provides a promising way 
to measure peculiar velocity at cosmological distance and 
makes itself a powerful probe of the dark universe. It has 
allowed the measurement of the structure growth rate in 
spectroscopic surveys such as 2dF d, |1] , SDSS l^, [lH , 
VVDS [11, WiggleZ [H, [l| and BOSS [H Such 
growth rate measurement is highly valuable in probing 
the nature of dark energy [17h24| and gravity |25l - l32| . 
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In particular, with both the expansion rate measurement 
from BAO |33l - [36| and structure growth rate measure- 
ment from RSD, spectroscopic redshift surveys are well 
suited to test consistency relations in general relativity 
and to discriminate between dark energy and modified 
gravity [s^- For these reasons, RSD has become one of 
the key science goals of the planned stage IV dark en- 
ergy projects such as the BigBOSS experiment [ssj ] and 
the Euclid cosmology mission [s^. It can also be used 
to measure the stochastic galaxy bias [40j and the galaxy 
(pairwise) velocity dispersion [41|, both are valuable for 
studying galaxy formation. 

These applications heavily relies on the RSD model- 
ing. However, modeling RSD to accuracy matching stage 
IV dark energy projects is very challenging, especially 
due to three sources of nonlinear and the associated non- 
Gaussianity. (1) The nonlinear mapping from real space 
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clustering to redshift space clustering. Due to this non- 
linearity, even the lowest order statistics in redshift space 
involves correlations between density and velocity fields 
to infinite order (e.g. '42'-'45]). (2) The nonlinear evolu- 
tion in the real space matter density and velocity fields. 
(3) The nonlinear and nonlocal galaxy-matter relation. 
The galaxy density bias is known to have non-negligible 
nonlinearity and stochasticity (e.g. ^M,)- The galaxy 
velocity bias may also be needed for precision modeling 
(e.g. [43] )■ To proceed, layers of approximations and 
simplifications have been made. 

Here we elaborate on some of these approxima- 
tions/simplifications with example of matter clustering 
in redshift space. One of the most commonly used RSD 
formulae connecting the isotropic real space power spec- 
trum Pss (fc) to the anisotropic redshift space power spec- 
trum Pis{k,u) is Pis{k,u) ~ Pss{k){l + fuYD''°^{ku). 
It is a phenomenological combination and extension of 
the linear Kaiser effect and the finger of God (FOG) ef- 
fect due to random motions. Here, / = din D/ din a and 
D = D{z) is the linear density growth factor at redshift 
z — 1/a — 1. Throughout this paper, the superscript 
"s" denotes the property in redshift space. Pgg depends 
on both k and u = k\\lk. Throughout the paper we 
adopt the z-axis as the line of sight, so /c|| = kz- Ap- 
proximations/simplifications made include (1) the par- 
allel plane approximation [4^, (2) linear evolution in 
the velocity-density relation and (3) no stochasticity be- 
tween the velocity and density field [l^. (4) It also 
neglects most high order correlations between the den- 
sity and velocity fields [H, Hil. (5) £)^°°(fcu) is a phe- 
nomenological description of the FOG effect, the overall 
damping caused by random motion. Both a Gaussian 
form D^'~"^{ku) = exp{—{kua.u/H)'^) and a Lorentz form 
L'™°(/cm) = 1/(1 -I- {kuap/Hf/2) are widely adopted. 
However, the physical meaning of the velocity dispersion 
ay and especially the pairwise velocity dispersion dp is a 
bit ambiguous. Furthermore, deviations from the above 
forms have been found in simulations [49j . 

Various approaches have been investigated to improve 
the RSD modeling. A far from complete list includes the 
Eulerian and Lagrangian perturbation theory [4^ . l47l [50l- 

f, the halo model IH, |57| - |60| . the streaming model 
[6ll . [63. t he recently proposed distribution function 
approach [4a - |45| and combinations between them. Fur- 
thermore, due to significant nonlinearities involved, RSD 
modeling often resorts to numerical simulations on cal- 
ibration and testing (e.g. [H, [H, [H IH-liil). Despite 
these efforts, RSD modeling has not yet achieved the ac- 
curacy required for precision RSD cosmology. For exam- 
ple, recent tests against N-body simulations find that the 
inferred / can be biased low by 0(10%) or more [63l - [68j . 
significantly larger than the 0(1%) statistical error in / 
for surveys like BigBOSS and Euclid. 

Since RSD is induced by peculiar velocity, it is of cru- 
cial importance to understand the peculiar velocity field. 
We find that, an appropriate velocity decomposition has 
the potential to simplify and improve the RSD model- 



ing. We also find that the same decomposition may 
enable robust reconstruction of 3D peculiar velocity in 
spectroscopic surveys. It has the advantage to render 
the otherwise troublesome RSD in velocity reconstruc- 
tion into valuable source of information. The current 
paper aims to lay out the methodology of the proposed 
velocity decomposition, RSD modeling and 3D velocity 
reconstruction. Extensive tests against simulations and 
mock catalogs are required to clarify /justify/quantify nu- 
merous technical issues. These numerical results will be 
presented in companion papers [l| . 

This paper is organized as follows. In fJTT]we decompose 
the peculiar velocity fields into three eigen-modes (v^ , vs 
and Vs) and discuss related statistics. Among them, 
is the velocity component completely correlated with the 
density distribution and is the one contains most cos- 
mological information. In i jllll we show that the three 
velocity components affect RSD in different ways and 
their impacts can be treated separately. We are then 
able to derive the exact RSD formula and, furthermore, 
make reasonable approximations for realistic calculation. 
Through this methodology, we explicitly identify a signif- 
icant systematical error in RSD cosmology. i jlVI proposes 
a method to reconstruct the 3D from redshift surveys 
and the appendix [D] proposes an alternative. For brevity, 
the above sections focus on the matter field. But these 
results can be extended to the galaxy field straightfor- 
wardly, as briefiy discussed in |jVl We further argue that 
the velocity reconstruction is insensitive to the galaxy 
bias. At the end of each sections, we list key statistics 
to be investigated in future works. We also prepare four 
appendices for technical issues. 

II. PECULIAR VELOCITY DECOMPOSITION 

Any vector field can be decomposed into a irrotational 
(gradient) part and a rotational (curl) part. Analogous 
to the electric and magnetic fields (and also the CMB 
polarization field and the cosmic shear field), we denote 
the first one with a subscript "E" and the later one with 
a subscript "B". Hence the peculiar velocity v can be 
decomposed as 

V(x) =V£;(x)+Vb(x) . (1) 

By definition, V x = and V • = 0. In Fourier 
space, we have V£;(k) = [v(k) • k]k and VB(k) = v(k) — 

V£(k). 

can be completely described by its divergence 
6'(x) = -V • v(x)/iJ = -V • vs(x)/if. To conveniently 
describe the velocity-density relation, we carry out a fur- 
ther decomposition, 

V£;(x) =V5(x)+Vs(x) . (2) 

Both V5 and V5 are irrotational (V x — and 
V X V5 = 0). We require that 9s = ■ vs/H is com- 
pletely correlated with the overdensity S. So we denote 
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this component with a subscript "5". 
we then have 

0s{k) ^ S{k)W{k) , 



In Fourier space, 



(3) 



where PF(k) is a deterministic function of k, to be deter- 
mined later immediately. 

On the other hand, we require that 6s = —V • vs/H 
is uncorrelated with S ((6's(x)<5(x + r)) — 0). 6s is the 
source of stochasticity in the 5-6 relation, so we denote 
this component with a subscript "S" . 

We define the power spectrum between field A and 
field B through {A{k)B{k )) ee (27r)3^3D(k + k')PAB(k). 
We often use the notation A\g = Pab / C^tt^) , 
which enters into the ensemble average (A(x)B(x)) — 
J A\B{k)dk/k. Through the relation {S{k')6{k)) = 
{d{k!)0s{k)) = {S{k')S{k))W(k), we obtain 



W{k) = W{k) = 



Psejk) 
Pssik) 



(4) 



Notice that, due to the isotropy of the universe, Pse{k) — 
Pse{k) and Pssik) = Pss{k). So W does not depend on 
the direction of k. 

Hence we prove that the velocity decomposition into 
V5, V5 and is mathematically unique, with no as- 
sumptions on the underlying density and velocity field. 
But the decomposition has strong physical motivation 
too. The three components are associated with different 
physical processes in the structure formation. We now 
proceed to their physical meanings. 



A. The field 

In the limit k fcNL, V5 is the only velocity compo- 
nent. Here ^nl is the nonlinear scale, defined through 
A^i-(fcNL) = 1- is a decay mode before shell cross- 
ing so it is negligible in the linear regime [6^ [T^I- On 
the other hand, in the limit k <^ ^nl, is completely 
correlated with the density field, with a deterministic re- 
lation 6 = fS. So V5 vanishes too. v^, being the most 
linear velocity component, is of the greatest interest to 
cosmology. 

Nevertheless, nonlinear evolution leaves non-negligible 
imprints in the field. A crucial point is that nonlinear 
evolution affects the density field and the velocity field 
in different or even opposite ways (e.g. [i^, HI])- The 
third order Eulerian perturbation shows that, when the 
effective power index Ucff ^ —1.9, nonlinear evolution 
actually suppresses the velocity growth, while enhances 



the overdensity growth (Fig. 12, [69|). Even for 



< 

— 1.9, 6 grows more slowly than the overdensity (Fig. 
12, [69|). Furthermore, nonlinear evolution generates a 
stochastic velocity component vg, further reducing 6s 
with respect to S. In the deeply nonlinear regime, after 
many orbit crossings, the velocity field eventually loses its 
correlation with the density field and we expect — 0. 
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FIG. 1: Solid lines are W{k,z) predicted by the third 
order Eulerian perturbation theory at redshift z = 
0.0, 0.5, 1.0, 1.5, 2.0 (bottom up). We adopt a flat ACDM cos- 
mology with rim = 0.26, fif, = 0.044, flA = 0.74, h = 0.71, 
as = 0.8 and Us — 1. Nonlinearity drives W down from 
unity. Later we will show that what inferred from redshift 
space distortion is fW. So the structure growth rate / can 
be biased low by 0(10%), consistent with recent findings 
(e.g. [g^ - Ig^ I. To demonstrate the vs component, we also 
plot yj Pee / Pss, normalized to unity at fc — >■ (dash hues). 

yPeeJPss = Pse/Pss x VI + V- V = Peses/PosOs quanti- 
fies the relative amplitude of vg with respect to the veloc- 
ity component vj. It also quantifies the velocity divergence- 
density stochasticity (Eq. [T]). In the limit fc — >■ 0, 77 — >■ and 
yPee/Pss — > Pse/Pss- The nonlinear evolution generates vs 
and causes the two sets of curves to deviate from each other. 
Since the third order Eulerian perturbation theory has lim- 
ited range of applicability, numerical results shown in this 
plot are mainly presented to demonstrate the major impacts 
of the nonlinear evolution, to drive W down from unity and 
to drive rj up from zero. Robust quantifications of W and 77 
will be presented in companion papers 1] . 



Hence we should use Eq. U] to describe the 6s-6 re- 
lation, instead of the linear relation ds = fS. For the 
convenience of highlighting the deviation from the linear 
relation, we define the normalized W through 



W{k) 



W{k) 



Psejk) 

Wik^O) f Pssik) 



(5) 



W calculated using the third order Eulerian perturba- 
tion theory is shown in Fig. [TJ In companion papers 
[ij we will numerically evaluate this key quantity using 
high resolution N-body simulations. In the limit fc — > 0, 
1 as expected. But nonlinear evolution soon drives 
W to deviate W < 1. Even at relatively high redshift 
z = 2 and pretty linear scale k = 0.1/i/Mpc, the devi- 
ation already reaches 0(1%). The derivation increases 
towards low redshift and exceeds 10% at z < 0.5. 
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To better understand this behavior, we can express W 

as 

Here, rse{k) is the cross correlation coefficient between 
S and 9. By definition, rse{k) < 1. Given the fact that 
velocity growths more slowly the density, Pgg < pPss- 
So we expect W < I'm nonlinear regime and expect that, 
when k ^ /cnl, fse — > and — > 0. 

The recognition of the velocity component has sev- 
eral important applications. (1) It lays out a promis- 
ing way to reconstruct the 3D peculiar velocity from 3D 
density distribution, which is accessible from galaxy spec- 
troscopic redshift surveys. W mimics a window function 
with a smoothing scale comparable to the nonlinear scale. 
It exerts on the density field and suppresses small scale 
inhomogeneities so that the smoothed density field pro- 
vides a honest estimation on the underlying velocity field 
v^. One can find a similar window function exerting on 
the redshift space density, which is directly observable. 
Hence the reconstruction of the stochastic 3D vector field 
is simplified as the reconstruction of a deterministic win- 
dow function with limited degrees of freedom. This is 
one of the most important applications of the velocity 
decomposition proposed in this paper. Later in i jlVI we 
will present more detailed investigation. (2) It simplifies 
the RSD modeling. Since — J> toward small scales, 
nonlinearity and non-Gaussianity in the field are sig- 
nificantly suppressed. (3) It identifies a severe system- 
atical error in RSD cosmology. Later in SJHTC] we will 
show that RSD is determined by fW instead of /. If the 
factor W is not included in the RSD cosmology, / will be 
biased low by ~ 2% at z = 2 and ^ 10% at z = 0, even 
if we restrict the analysis to fc < 0.1/i/Mpc. To our best 
knowledge, this is the first time such systematical error is 
diagnosized explicitly. This source of systematical error 
could explain recent findings of 0(10%) underestimation 
in / inferred from simulated RSD data [gs'Jgs']. We will 
present more discussions on its cosmological implications 
in TOEI 

B. The vs field 

To the opposite of v^, V5 vanishes in the linear regime 
and begins to grow due to the nonlinear evolution. For 
this reason, it lacks large scale power and hence its cor- 
relation length is smaller than that of . For the same 
reason, it is intrinsically non-Gaussian. Being uncorre- 
lated to the density field, V5 induce stochasticities in the 
velocity-density relation. We have the relation 

rsBik) = , = , , (7) 

^P55{k)Pee{k) ^T+W) 

where ■q{k) = Pg^gg{k)/ Pggg/j{k). When V5 overwhelms 
over Vs (ri ^ 1), rsg and the velocity field loses its 



correlation with the density field. 

Based on the third order perturbation calculation, we 
find that PggSs already reaches ~ 1% of Pg^es (namely 
ij ~ 1%) at = 0.1/1/ Mpc and z ^ 1, as can be in- 
ferred from Fig. [TJ As expected, the situation is less 
severe at higher redshifts. But even at z = 2, 77 ~ 1% at 
k = 0.2/i/Mpc. Hence in general, V5 is a non-negligible 
velocity component even at scales which are often con- 
sidered as linear. Nevertheless, it is still subdominant to 
Vs at scale k ^ 1/i/Mpc, as one can infer from Fig. 1 
of p^ . in combination with our Eq. [T] Our numerical 
evalulations obtain similiar results and we will present 
details in companion papers [l| . 

What cosmological information does it contain? One 
particularly interesting piece is likely the nature of grav- 
ity. For modified gravity models to pass the solar sys- 
tem tests while driving the late time cosmic acceleration, 
gravity must behave upon the environment (e.g. review 
articles [7ll,[73|). Such environment dependence often be- 
comes prominent in the nonlinear regime and can signif- 
icantly affects the velocity divergence (e.g. [s^l). Being 
arising from nonlinear evolution, vg would be sensitive 
to this smoking gun of modified gravity, making it at- 
tractive for testing gravity. 

In j jllll we will show that the vs induced RSD differs 
from the induced RSD (i jIII[) . We argue that these 
differences can be used to separate the different veloc- 
ity components and made statistical measurement of V5 
possible. 



C. The Vfl field 

vb decays as long as the single fluid approximation 
for the dark matter distribution holds (69'|. vb grows 
only when the nonlinearity is sufficiently large that shell 
crossing happens. In the deeply nonlinear regime, we 
may even expect equi-partition in the velocity distribu- 
tion and Vb can dominate over ve- So we expect that its 
power concentrates at even smaller scales than vg. In- 
deed, numerically studies [t^ show that the power of this 
velocity component concentrates at small scales, with 
r.m.s. much smaller than that of v^. Our numerical 
evalulations obtain similiar results jlj]. Later in mil Al we 
will show that these behaviors significantly simplify the 
modeling of induced RSD. For in depth study of v^, 
we refer the readers to (70j and references therein. In [ij , 
we will present the statistics of vs in the context of RSD. 



D. Statistical description of the three velocity 
components 

We are not at a position to calculate the statistics of 
these velocity components. Instead, we present some 
general discussions here and postpone any quantitative 
calculations in future works. First, since the three ve- 
locity components have different origins, the halo model 
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scenario [41|, |73| - |75| provides an unified way to describe 
the three components, with a number of extensions. (1) 
Since deviations from the linear theory prediction of the 
halo bulk motion have been diagnosized [76| and since it 
fails to describe the emergence of vs, a natural tool to 
describe the halo bulk motion is the third order pertur- 
bation. (2) Peculiar velocity has been treated as the sum 
of the halo bulk velocity and the random velocity inside 
of the virialized halos [73] . It may be extended to include 
the more complicated motion around the halo outskirt, 
which may be a significant source of [t^I . Neverthe- 
less, treatment of is trickier and we refer the readers 
to [20] for detailed discussions. 

One important statistics relevant to RSD is the ve- 
locity correlation function. Due to symmetry considera- 
tions, it can be decomposed as [4j 

K(xi)i,,(x2)) = ^|JAr)S^, + [^11 (r) - V^W] ^ . (8) 

Here r = Xi — X2 is the pair seperation vector. is the 
correlation function when both Vi and Vj are along r. 
is the one when both velocities are perpendicular to r. 

i/'ll and ip± are not independent. For a potential flow 
{vs and V5), we have the textbook result jj] 



d{ril;j_{r)) 



dr 



2 sin(fcr) dk 



(9) 
(10) 



Here we have defined the velocity power spectrum 
through (v(k) •v(k')} = (27r)3<53c(k k')P™(fc)- Its 
covariance is defined as = P^^k^ /2ir^ . 

For Vs, we no longer have this relation. Using the fact 
that Vs can be expressed as the vorticity of a vector A 
(vs = V X A), we derive the following relations. 



'0±(r) 
V'll W 



1 rfV'ii(r) 



dk 1 
T3 



sin kr cos kr 



{krY [krY 



(11) 
(12) 



E. To do list 

In companion papers [l|, we will use N-body simula- 
tions to numerically evaluate statistics of these velocity 
fields. An incomplete list includes 

• The power spectra P^^^_5(fc,z), Pygv^{k,z) and 
Pygyg{k, z). Their scale and redshift dependences 
are helpful to understand the impact of nonlinear 
evolution. 

• The correlation function 'il}y^v^{r, z), 'ipv^vc:{r, z) 
and 'ipvgvg{r, z), of the "||" mode and the "±" 
mode. These statistics quantify the correlation 
length of the three velocity components and tell 
us at which separation we can treat the velocities 
at two positions as independent. 



The PDFs of v^, vs and v^. The PDFs determine 
the overall damping (FOG) to the redshift power 
spectrum. We will also calculate the cumulants to 
quantify non-Gaussianity of hese velocity compo- 
nents. 

W{k, z). As we have addressed, W is of crucial im- 
portance in the 3D velocity reconstruction and in 
the RSD cosmology. Through N-body simulations, 
we will robustly measure its dependence on k and 
z. Furthermore, we want to quantify the W{k, z)- 
A'gg{k,z) relation to better understand its depen- 
dence on the nonlinearity. Especially, we want to 
know if it can be well approximated by a simple 
function with only a few parameters. 



III. MODELING THE REDSHIFT SPACE 
DISTORTION 

Now we proceed to the RSD modeling with the aid 
of the proposed velocity decomposion. Following 42], 
we utilize the matter conservation ((1 -I- (5'*(x'*))(i^x'* = 
{1 + 6{x))d'^x) to derive 



S%k) = y" [1 + ^W] exp (i'^^ exp(zk • x 



)d^x (13) 



In the above expression, we have neglected a Dirac fun- 
tion which only shows up when k = 0, irrelevant to 
our calculation. This equation adopts the plane par- 
allel approximation and adopts the line of sight as the 
z-axis. It also assumes no multiple streaming. Exten- 
tion to this more complicated situation can be done by 
the phase space distribution function approach [43-45]. 
Phenomenologically speaking, multiple streaming can be 
described by an overall damping to the redshift space 
power spectrum. Hence we do not expect major changes 
to the results presented in this paper. 



The redshift space power spectrum is given by 42 1 



<5<5 



1 + (5i)(l + ^2) exp A) exp(-ik ■ r)d^r{U) 



where Si = 6{Ki) and X2 = Xi -|-r. The ensemble average 
(• • •) is averaged over all xi with fixed r. For brevity, 
we have denoted A = ik^iviz — V2z)/H. We also denote 
Aq = ikz{viz,a — V2z.a)/H whcrc a = S,S, B. Due to the 
axial symmetry along the line of sight, P^^ only depends 
on k and u. Hereafter we often write it as Pgsik, u). 

An immediate result is that 5^ (k_L ,kz — 0) = 
Sik±,kz = 0) and Pisik±,kz = 0) = Pss{k±). A sim- 
ilar relation holds for the bispectrum, i?|(ki, k2, ka) = 
_B3(ki, k2, ks) when ki^z = ^2,2 — ks^z = 0. Later we will 
show that these relations are very useful in RSD model- 
ing. Notice that no such relations exist in real space. For 
example, £.^{r±,rz = 0) ^ i{r±)- 
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A. The role of vb 

Now we will apply the velocity decomposition to Eq. 
1141 Since only is correlated with the overdensity, Eq. 
[TH can be reduced to 



P!sik,u) = j {{il + 5,){l + S2)e^pXs) (15) 

X Ds{k,,r)DB{k,,r)} exp{-ik ■ r)d^r . 

The two function Ds{kz,r) and Dsikz, r) completely de- 
scribe the redshift distortion caused by V5 and respec- 
tively. 



Ds{kz,r) = ( exp U 



■ kz{viz,R - V2z,r] 



H 



DB{kz,v) = ( exp 



kzjvizM ~ V2z,b) 

H 



(16) 



(17) 



Due to the symmetry between kz and —kz^ Ds.sikz, r) = 
DsM\kzlv). 

We have Ds < 1 and Db < 1, so both V5 and 
suppress the clustering presented by terms inside of the 
bracket of Eq. [131 However, V5 and vb can have their 
own clusterings and hence can in principle increase the 
overall clustering in redshift space. We can also define a 
corresponding function for v^, 



Ds{kz,r) = (exp 



. kz{viz,5 - V2z,5) 

H 



(18) 



But since is correlated with the density field, its role 
in RSD is too complicated to be described by a single 
function Dg- 

In general, these functions depend on the pair separa- 
tion r, due to correlations in the corresponding velocity 
fields. To describe this effect, we define 



l + e5,S,B(r, fcz) 



Ds..s,B{kz,'r) 
Ds.s,B{kz,r 00) 



(19) 



Here, the limit r ^ 00 corresponds to the limit of no 
velocity correlation. In this limit, peculiar velocity only 
causes damping (FOG). Hence we denote D at this limit 
with a superscript "FOG" , 



exp 



kzV, 



z Vz,5,S,B 

n 



As discussed in []TT1 arises mostly in the deeply non- 
linear regime, so we expect a correlation length shorter 
than scales of interest for RSD cosmology. Hence we can 
neglect correlation in the field. Namely, we can make 
the approximation 







DB{kz 



nFOG 

^B 



{kz 



(21) 



Since D^'^ is independent of r, Eq. [15] is simplified as 



P!5ik,u) = |y {{l + 6i){l + S2)expXs) (22) 

I?s(fc.,r)exp(-ik-r)d3r}i?™G(fc,) . 

The damping function D^'~"^ is related to its velocity 
PDF through 



/OO 
exp 
-00 



.kzVz,B 1 „ . . , 
« -J^ ) r'B[Vz^B)dVz,B 



kzVz,B In/ N J 

COS I — ) FB{Vz^B)dVz,B 



Through the cumulant expansion theorem, we can ex- 
press Z?^*-*^ in cumulants of Vz.b- The cumulant expan- 
sion theorem states that 



exp I 



kzV 
H 



= exp 



I 2 2 \ / °° 

^ Vj>2 



exp 



2i 



2j 



12' 



2 

360^ 



Here, v denotes Vz^b- (• • •)c is the corresponding cumu- 
lant of the quantity inside of the bracket. For a Gaussian 
velocity distribution, we have (w-')c — for j > 3 and 
recover the Gaussian FOG. For non-Gaussian distribu- 
tion, higher order terms show up. Since (f^-'^^)c = 0, 
only even cumulants contribute. Kj = {v^)c/(tI is the 
reduced cumulant and x = {kzC^/H)^. 

Non-Gaussianity in could be significant and we may 
expect that high order cumulants must be included to 
robustly model D^'-"^. Fortunately, the reality can be 
much simpler due to the fact that vs has small r.m.s. 
Its cr„ is more than a factor of 10 smaller than that of 
V£; [2^. This means that < 10 km/s even at z = 0. 
Hence for the scales of interest (e.g. k < l/i/Mpc), we 
have X ^ 10~^. Then the induced FOG can be well 
described by the following Gaussian form, 



riFOG 
^B 



{kz 



exp 



(24) 



(20) 

is hard to calculate from first principle, so it shall be 



treated as a free parameter to be fitted by the data. 



B. The role of vg 

The situation for vg is more complicated. Since a sig- 
nificant fraction of V5 comes from bulk motion, V5 can 
be still correlated over 0(10) Mpc separation. So we are 
no longer able to make the approximation £5 = 0. In- 
stead, we have, with the aid of the cumulant expansion 
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theorem, 



1 + es(r, /e^) = exp 



kl{viz^sV2z,s) 



(25) 



Xe['='<3t-L,s''L,s-2(fL,s»'2.,s + l«2))J12if* + ---] 

Due to the symmetry v o —v, odd terms in v do not show 
up in the above equation. Plugging it into Eq. 1221 we 
find that the leading order contribution to is Pg^g^u'^ 



(Eq. which enhances P|^. Notice that, since vg is 
uncorrelated with the density field, it does not contribute 
a term as does. Higher order contributions come 
from both fourth and higher order velocity correlations 
(Eq. [^5)1 and from the convolution of eg with ((l-|-(5i)(l + 
52)exp\s) (Eq. [T5]) . These corrections are of the order 
5'^ or higher. We group all these contributions into a 
single function Cs{k,u) (Eq. I^S)) . Later we argue that 
we may be able to set Cs = in the RSD modeling. 

On the other hand, V5 causes an overall damping, char- 
acterized by Df'-"^, 



DlOG^ku) = P exp [i^^) Ps{v..s)dv.,i 



cos 



V H 



Psivz,s)dv^^s 



Since vg arises from nonlinear evolution, it is intrinsi- 
cally non-Gaussian. We will check if D^'^'^ can be well 
described by the first several cumulants, with the aid of 
the cumulant expansion theorem (Eq. 



C. The role of Vi 

Due to coupling between the velocity and the density 
fields, redshift distortion induced by is complicated, 
although analytical expression in the Gaussian limit ex- 
ists [42| . In the appendix, we derive the non-Gaussian 
corrections (Eq. IA4p . In Fourier space, we have 



55 



{k,u) = {Pss{k){l + fW{kyf 



u'^Pf. 



OsOs 



(k) 



-CNG{k,u) + CG{k,u) + Cs{k,u)} 



xD 



FOG/ 



(fcu)Z?|°^(fcu)i?^^^(fcM) . (26) 



FOG, 



Here, Cg and Cng are the Fourier transforms of 
CG(r, fc;,) (Eq. \M\l and C7VG(r,fcz) (Eq. EZl), respec- 
tively. Cg contains all high order corrections if the den- 
sity and velocity fields are Gaussian. Cng is the extra 
correction caused by non-Gaussianity. Both terms are 
consequeces of the vs-S correlation. 

vs also causes a FOG effect, described by Dj'~^'^. 
is largely Gaussian due to the suppression of < 1 in 
the nonlinear and non-Gaussian regime. However, due to 
the large amplitude of (t„^ , it is unclear whether we can 



neglect the K4 and/or Kq corrections. Here cr„^ is the 
the one dimension velocity dispersion of v^. 



2 _ 



- 3 [ KvAk)— 



(27) 



1 1 ^^Uk)w^k)'± 



The terms CnGi C'g and Cs are the sums of infinite 
power series of 5. To carry out realistic calculation, we 
need to truncate them in a reasonable way. The first line 
terms in the right hand side of Eq. [55] exhaust contri- 
butions of the order 5'^. It can be rewritten in a more 
famihar form Pgg + 2v?Psg + u'^Pge (e.g. [12]) and con- 
firmes previous results. The second line terms are higher 
order in S. The leading order term of Cng is cx 5^, while 
those of Cg and Cs are cx (5*. In this sense, for Eq. [26] 
to be complete at the order of , we can set Cg = and 

C5 = 0. 

But in term of the linear density 6l, the situation is 
different. (1) PssOs vanishes in the linear perturbation 
theory, but emerges in the second or higher order Eule- 
rian perturbation theory. So Pg^Os itself is of the order 
Sj^ and Cs is of the order (5£ . (2) In the Eulerian pertur- 
bation theory, leading order terms of Cg and Cng are of 
the order Sj^. So up to the order (5|^, we can set Cs = 0. 

Hence no matter in power series of S or in power series 
of Sl, we can set Cs = when calculating the leading 
order corrections to the usual RSD formula. But whether 
or not we shall set Cg = is an issue for numerical 
examination. 

The inclusion of Cg and Cng may appear to introduce 
more difficulties and uncertainties to the RSD modeling. 
However, this is not the case. We will show that, the in- 
clusion of Cg and leading order terms in Cng does not 
introduce extra degrees of freedom in the RSD modeling. 
So the inclusion of these terms is capable of reducing sys- 
tematical errors, without degrading cosmological param- 
eter constraints. This is definitely a desirable property, 
made possible by the proposed velocity decomposition. 



1. The Cg correction 

C'g is an analytical (but nonlinear) function of the two 
point density and velocity (v^) correlation functions (Eq. 
IA5|) . Since Pgs is directly measurable from the fc^ = 
Fourier modes in redshift surveys, Cg up to any order 
can be calculated strictly without introduce any extra 
parameters. 

Cg does not contain terms of odd order in 6. So we can 
spht Cg as CG{k, u) — J2j CG,2j{k, u) where j = 2, 3 • • •. 
In reality, we may only need the leading order 



CG,4:ik,u) 



(2^ 
G(k,ki, 



Pssiki)Pssik2) 

k2). 



k^,k2,k''M^W,. 



(28) 
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Here, = k + ki and Wi = W{h). The kernel G is 



G 



k2zW2/kl 



kl 



2A^2 ^2 



(29) 



With an observable and W a function to be fitted 
anyway, calculating Cg requires no extra free parameters 
and hence does not induce new model uncertainties. 



2. The Cng correction 

CNoik^u) ~ J2j>3^NG,j{k,u) contains connected 

part of j-th order correlation ((v;^), {Sv-^'^) and (5(5w^~^)) 
to J — ?• cxD (Eq. IA7[) . These terms are the consequence 
of nonlinear and non-Gaussian evolution, so we expect 
them to become non-negligible only in nonlinear regime 
and may become dominant in deeply nonlinear regime. 
However, since — > in deeply nonlinear regime, their 
contributions to Pgg are significantly suppressed. For the 
same reason, Cngj+i is suppressed by a factor ~ W, 
with respect to Cng,]- Hence we propose to keep only 
the J = 3 term and neglect higher order corrections. 

Fourier transforming CArG_3(r, k^) (Eq. IA7I) . we obtain 



D. A new formula on the redshift space power 
spectrum 

We then propose the following formula for the redshift 
space power spectrum. 



P'{k,u) 2, I^Pss{k)(l + fW{k)u^y 



(31) 



-u^Pesdsik) +CG{k,u) + CNGAf^,u)} 



In Eq. EH W{k) = fW{k) and PesOsik) are the cos- 
mological information we seek for. We address one more 
time that, Cg and Cisig,3 oltr not free functions. They 
are uniquely determined by W . This valuable property is 
achieved by the proposed velocity decomposition and, in 
particular, by the deterministic relation 9s — 5W . 

There are only very limited degrees of freedom in the 
FOG terms. Through Eq. [H Dj^'^{ku) and Df^'^{ku) 
may be well described by a^^, Cug, K^^y^ and/or 

Ke,vsi ^6,Ds- Among them, a^^ is determined by W 
(Eq. [77)1 and a-u^ is determined by PesOs- So in principle 
they are not free parameters. Furthermore, since cr^^ ^ 
lOOcr^^ [tJ, we may be able to set ay^ ~ 0, depending 
on the desired accuracy. 



CjVG,3(A,«) = 2 



^B3(ki,k2,k)-^VF(fci) 



{2n) 



§Wik) + '-pWik2) 



1 



.(30) 



Here, k2 = — ki — k. i33(ki, k2, ks) is the real space mat- 
ter bispectrum. We notice that the ensemble average of 
B^Cki, k2, ks) is directly available from the same redshift 
surveys used for RSD measurement. Since ^ ■ kj = 0, 
kl, k2 and ka lie in the same plane. Due to the isotropy 
of the universe, i33(ki, k2, ka) does not depend on the 
inclination of the plane. So its value is equal to the 
case where all k lie in the plane perpendicular to the 
line of sight (namely the x-y plane). Redshift distor- 
tion does not alter the kz = Fourier mode. Namely, 
(5''(k^, kz = 0) = 6(k±, kz = 0). This means that can 
be directly measured from the observed Fourier mode 
with ki^z = 0. Hence including Cng,3 in the calculation 
does not introduce extra fitting parameters and hence 
does not weaken the cosmological constraints. 

The same trick does not apply to CNG,j>i- They in- 
volve 4-th or higher order correlation ((5(ki), • • • , (5(kj)), 
with kl + • • • -I- kj =0. In general, ki(i = 1, 2, 3, 4 • • • , j) 
do not lie in the same plane. We are no longer able to 
infer their values from the kz — i) modes. However, due 
to extra suppression caused by < 1, we do not expect 
these terms to be important. Nevertheless, the accuracy 
of neglecting these higher order terms must be quantified 
through N-body simulations. 



1. Comparing to existing models 

Here we compare the proposed RSD formula (Eq. I3ip 
with a few existing RSD formulae. We choose these for- 
mulae because they can be compared with ours relatively 
straightforwardly. So this comparison is by no means 
complete. For brevity, we focus on the matter power 
spectrum for which we do not need to worry about the 
galaxy bias, especially the nonlinear/non-deterministic 
bias. (1) The Kaiser formula plus the FOG effect (here- 
after KF), 



Pls{k,u) ~ Ps5{k){l + fu'YD''''-{ku) 



(32) 



This is perhaps the most commonly used RSD formula. 
(2) The Desjacques & Sheth 2010 formula (hereafter 
SDIO, [13), 



Pls{k,u) = Pssim + fu^f 

xexp{-klal/H^)V^i,{kz) 



(33) 



Here, I/vir is the damping caused by random motions of 
virialized particles in halos. ay is the velocity dispersion 
of bulk motion. It improves over the KF formula by 
correctly recognizing the two velocity components (bulk 
motion and random motion). (3) The Scoccimarro 2004 
formula (hereafter S04, [42]), 



PMk,u) = {Pss{k) + 2u^Psgik) 
X eM-k^^v/H^) ■ 



'Peeik)) 



(34) 
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Here, can be treated as the one calculated from the 
perturbation theory or as a free parameter. It improves 
over KF by taking into account differences in the den- 
sity and velocity nonlinear evolution. As a consequence, 
it does not assume a deterministic relation between den- 
sity and velocity. (4) The Taruya et al. 2010 formula 
(hereafter TIO, [111), 



of the ACDM cosmology and confirmed its validity [78|. 
Given the existence oiW^l, the expectation value of 
Eg should be corrected by a factor 1/W, 



(35) 



+ A{k,u) + B{k,u))exp{~k^u'^(Tl/H^ 

The extra term A involves integral over the bispectrum. 
The term B involves convolution of two power spectra. 

Comparing to Eq. 1311 we recognize a number of cor- 
rections to existing formulae. Here we briefly summarize 
the most significant ones. (1) Both KF and DSIO fail to 
capture the W correction on /. So / based on KF and 
DSIO is underestimated by a factor W. (2) S04 and TIO 
improve over KF and DSIO by correctly capturing the 
W correction, although S04 and TIG do not make this 
correction explicitly. TIG further improves over S04 by 
including next order corrections (A and B). In the ap- 
pendix |B1 we prove that A — Cng.3 in the limit vs — ^ 0. 
However, TIO adopted the approximation = (equiv- 
alently {V1V2) = 0). so it fails to capture some terms of 
the same order as A and B, as we do (Eq. IA6p . More 
importantly, calculating A and B requires heavy model- 
ing or simulation calibration. Our Eq. [31] avoids these 
uncertainties. (3) The FOG effect is caused by three dis- 
tinctly different velocity components, so it may not follow 
a simple function form. For this point, the closest match 
to our formula is DSIO. These complexities can bias the 
interpretation of the inferred ay from RSD. 



E. Implications on the RSD cosmology 

Our analysis above shows that what one can infer 
from RSD is the combination fW{k) instead of /. Since 
W{k) 0.9 at fc = 0.1/i/Mpc and z — 0, this can cause 
10% underestimation in /. We expect that it is at least 
a significant source causing the recently found underesti- 
mation in / [63l - l68| . This systematical error overwhelms 
the 1% level statistical accuracy in / by stage IV dark 
energy projects. Our work provides a promising way to 
overcome this severe systematical error. To do so we 
need to robustly model W, with the help of analytical 
and numerical tools. 

The W correction affects many applications of redshift 
distortion. For example, ^25i] proposed an Eq estimator 
combining weak lensing and redshift distortion to test 
general relativity at cosmological scale. This estimator 
is insensitive to galaxy bias and is also less affected by 
initial fluctuations. [28| made the first Eq measurement 
and found Eq = 0.39 ± 0.06 at an effective redshift 0.32 
and O(10)Mpc//i scale. This measurement confirms gen- 
eral relativity within 20% observational uncertainties. 
This result has been used to perform consistency tests 



Eg = 



Gn f 



1 



W{k) 



(36) 



Future experiments such as BigBOSS+LSST, Big- 
BOSS-|-Planck CMB lensing, Euclid or other combina- 
tions of spectroscopic surveys and imaging surveys are 
capable of measuring Eq to 1% level statistical accuracy 
[25| . For these measurements. The 1/W correction is 
needed in order to correctly interpret the measured Eq- 



F. To do list 

Most statistics discussed in this section are too compli- 
cated to evaluate analytically. In companion papers [Ij], 
we will use N-body simulations to numerically evaluate 
these statistics. An incomplete investigation list includes 

• The damping function Dj'~"^{ku), D^'~"~^{ku) and 
D^'-"^{ku). In particular, we will investigate the 
usage of the cumulant expansion theorem (Eq. [23| . 
We will check if including K4 and/or Kq describes 
Dg'~"^{ku) and D^'~"^{ku) accurately at scales of 
interest. We will also check the accuracy of Eq. [511 

• es, es a-iid es- In particular, we will quantify the 
accuracy of Eq. [2T]&[25l 

• The accuracy of the proposed RSD formula (Eq. 
[3T|) . We will measure Cg.j and Cngj from simu- 
lations. We will also quantify the relative contri- 
bution of the leading terms, namely Co, 4 to Cg, 
Cng,3 to CjvG and PegBs "^ith respect to Cs- We 
will seek the possibility to further improve the RSD 
modeling, if needed. 



IV. RECONSTRUCTION OF THE 3D 
VELOCITY FIELD IN REDSHIFT SURVEYS 

So far we focus on inferring the statistical average of 
the velocity field including W oc Psg, Pvsvs^ '^vsi <^vs 
and . Can we go a step further to construct the full 
3D velocity field? The gain would be huge since the 3D 
velocity field contains much more information and has 
much more applications. One of such applications is to 
search for the missing baryons through the diffuse kinetic 
Sunyaev Zel'dovich (kSZ) effect [t^ caused by the inter- 
galactic medium. However, although the ACT collabora- 
tion has reported the first detection of the kSZ effect of 
galaxy clusters [131, the diffuse kSZ is still elusive even 
to the state of art surveys such as ACT [13, HI and SPT 
[s^ . [s^ l. due to various overwhelming contaminations. 
The measurement can be revolutionized by the kinetic 
SZ tomography [H, HH] ■ The key is to reconstruct the 
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galaxy velocity field. One can then construct the galaxy 
momentum field. It can then be correlated with CMB to 
measure the kSZ effect while automatically eliminating 
all contaminations of scalar type. A combination of Big- 
BOSS and Planck is promising to measure kSZ to better 
than lOcr at a number of redshift bins (85l |. 

How to reconstruct the 3D velocity field? If we can, 
which velocity component can be reconstructed? There 
are numerous works over a long history. We are not at 
a stage to overview these works. Rather, we outline our 
approach, based upon the proposed velocity decomposi- 
tion. 

Since observationally we only have the density field, 
which is a scalar field, the information budget does not 
allow us to reconstruct all the three v^, vs and vb vec- 
tor fields. However, since is completely correlated with 
the density field, it is promising to reconstruct from 
the observed density field. Eq. [3] guides us to propose a 
linear estimator for the 3D peculiar velocity reconstruc- 
tion through the 3D density field. Since we observe the 
redshift space overdensity instead, this linear estima- 
tor should operate on 5^(k), 



(37) 



We use the superscript "s" to denote properties in red- 
shift space. W corresponds to an anisotropic window 
function operating on the anistropic density field in red- 
shift space. Now the reconstruction of the stochastic 
3D velocity filed is simplified into th reconstruction of 
a deterministic function . Due to the axial symmetry 
along the light of sight, W^^(k) = W'{k,u). So it only 
has 2D degrees of freedom. These degrees of freedom are 
further limited by the asymptotic behaviors of u) 
at fcu — >■ and at fc ^ oo. 

We want the peculiar velocity estimator (Eq. [37)) to 
have no multiplicative error. This requires 



Wik) W{k) Pss{k) 



a(k) 



'(k)V^/.-(k) 



(38) 



W differs from the real space one by a direction depen- 
dent factor 1/C^ (k) . This factors arises because the red- 
shift space density 5^ is not completely correlated with 
the real space density 5. To better understand this point 
and to derive Eq. |38l we carry out a decomposition of (5^ 
into two parts, 

5^ = 51 + 5%. (39) 

We require 5^ to be completely correlated with the un- 
derlying overdensity 5 and hence to 9s. This is the part 
that we can use to recover the peculiar velocity. For this 
reason, we label this part with a subscript 5g is 

uncorrelated with i5 and 9s- It causes the stochasticity in 
5-5'^ relation and contaminates the velocity reconstruc- 
tion. So it is labelled with a subscript "S". Since 5y is 
completely correlated with 5, 



with the deterministic function to be determined. 
Through the relation ((5(k)5''(k')) = (5(k)(5^(k')) = 
((5(k)(5(k'))a(k), we obtain 



prs 

a(k)= 



(k) 



Pssik) 



'(k)i 



' Pisik,u) 
Pssik) 



(41) 



Here, P^^ is the cross power spectrum between S and (5*. 
rJI is the corresponding cross correlation coefficient. It 
has the asymptotic behavior r""* — 1 when fcu — > and 
j,rs _^ Q ^]-^gj^ /cu — j> cxD. Notice that r^'^ depends on both 
k and u. So do Cy and W . Hereafter we will write them 
as r^'^ikju), W^{k,u) and Cy{k,u) to highlight these de- 
pendences. 

Both Pss{k) and Pgg{k,u) are observables. To evalu- 
ate W'^ requires just one extra input, r'^'^ik, u). r*"** has a 
well defined asymptotic behavior r'"'* — > 1 when fcu — ?> 
and r^'^ — >■ when \ku\ oo. In the appendix [Cl we 
show that r'"'* is uniquely fixed by W and other quantities 
which can be inferred from Pssik, u). Hence the observed 
RSD allows us to figure out and then carry out the 
velocity reconstruction. In this sense, RSD is a source 
of information crucial for the velocity reconstruction, in- 
stead of source of noise as in many other apprroaches of 
velocity reconstruction. 

In the appendix [Dl we propose another approach to 
reconstruct the 3D velocity field. It shows more clearly 
the crucial role of RSD in the velocity reconstruction. 



A. Errors int the reconstruction and remodies 

A correct (fc, u) allows to reconstruct free of mul- 
tiplicative error. Unfortunately, additive errors persist, 
due to the stochastic component 6g. The reconstructed 
velocity divergence is 



9si^) 
Osi^) 



9sik)+9U^) 
S'sik)W'ik,u) 



(42) 



The additive error W5(k) — iH9gi'k)'k/k^ induces an ad- 
ditive error to the auto power spectrum directly mea- 
sured from 9s. So cosmological applications of the di- 
rectly measured auto power spectrum can be very lim- 
ited, if any. 

However, a good thing is that its cross power spectrum 
with the density distribution is still unbiased. 



Pse = Pse - Psi 



(43) 



For the same reason, the additive error does not bias the 
kSZ tomography [sj, [s^ and hence does not bias the 
effort to search for missing baryons. 

Since 9s is completely correlated with the matter den- 
sity, we can measure the auto power spectrum unbiasedly, 
through the relation 



5:(k) = 5(k)a(k) 



(40) 



p2 

Pss 



(44) 
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B. To do list 

Through N-body simulations, a number of key issues 
will be investigated in future works, 

• r^''{k,u), W'^{k,u) and Cv{k,u). We are then able 
to quantify the additive error in the reconstructed 
velocity. Furthermore, we want to understand the 
origin of (5|, which causes the stochasticity in the 
S-S'' relation and degrades the reconstruction per- 
formance. 

• The accuracy of Eq. IC6I to model r*"** . 

• The accuracy of the velocity reconstruction. It is 
determined by the accuracy of the inferred W from 
the observed RSD and the accuracy of the modeled 
r""*. Measurement errors in the galaxy distribution 
such as shot noise further complicates the recon- 
struction. We will take them into account for more 
realistic quantification of the reconstruction perfor- 
mance. 

• Further investigation of the velocity reconstruction 
approach proposed in the appendix |DJ 



V. FROM MATTER DISTRIBUTION TO 
GALAXY DISTRIBUTION 

In reality, we have 3D galaxy distribution instead of 
3D matter distribution. The above methods to model 
the matter redshift space distortion and to carry out ve- 
locity reconstruction can be extended to the 3D galaxy 
distribution straightforwardly and robustly. To do so, we 
need to decompose the irrotional galaxy velocity v^; into 
two eigen-modes such that one part is completely cor- 
related with the galaxy overdensity and the other part 
is completely uncorrelated. We can then replace quan- 
tities of the matter field in previous sections to that of 
the corresponding galaxy field. For example, W should 
be interpretted as W — Pge/Pgg where the subscript 
"g" denotes properties of galaxies. Now W = PW and 
/3 = //&g(fc— i>0). In fact, we can replace the subscript 
"(5" in previous sections with the subscript "g" to obtain 
results applicable to galaxies. 

Nevertheless, there is an important exception. The 
W in the Eq estimator (Eq. should be still inter- 
preted as that of the matter field instead of the galaxy 
field. By the construction, Eq oc Pgs/PgS j!!]. Here, S 
is the matter overdensity. We can decompose the galaxy 
velocity with respect to the matter field, instead of the 
galaxy number density field. Under reasonable assump- 
tions of PgSs = |86| and no galaxy velocity bias at 
relevant scales, Pgs/PgO = Pgs/{WPgs) oc l/ifW). This 
property is desirable since the correction term W can be 
robustly quantified through N-body simulations, free of 
uncertainties in modeling galaxy formation. 



This straightforward extension to the galaxy field fur- 
ther demonstrates the power of our approach. For ex- 
ample, it significantly alleviates the problem of galaxy 
bias in the velocity reconstruction. Even better, in the 
limit of a deterministic bias {5g — bgS), completely over- 
come it. Now W = Pge/Pgg oc l/hg. The reconstructed 
3D velocity is oc W5g oc 6°, independent of hg. This 
is also the case for the inferred velocity power spectrum 
Pe,e, = PgejPgg- 

In future works we will redo the analysis of the mat- 
ter density field for the halo number density. We will 
then proceed to mock catalogs of galaxies. Eventually 
we plan to develop efficient codes applicable to real data 
of spectroscopic redshift surveys. 



VI. DISCUSSIONS AND SUMMARY 

We have laid out the methodology to carry out 3D 
velocity reconstruction from 3D matter and galaxy dis- 
tribution. The method is based upon a velocity decom- 
position into three eigen-modes with physical motivation 
and of mathematical uniqueness. The same decomposi- 
tion also helps the RSD modeling and we have derived a 
new RSD formula. Through it we find that the inferred 
structure growth rate based upon some simplified ver- 
sions of RSD modeling can be severely underestimated. 
In a series of companion papers T] we will analyze N- 
body simulations to measure statistics of the three ve- 
locity eigen-modes, to test the accuracy of the proposed 
RSD formula and to quantify the performance of the pro- 
posed velocity reconstruction. 
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Appendix A: Calculating the induced RSD 

Throughout this section we only deal with the ve- 
locity component v^-. Hence for brevity we neglect 
the subscript "(5" and denote viz^s — ui, V2z.s — V2, 
Ai^2 = ikzVi 2/H and A5 = Ai — A2. We will adopt two 
tricks in [i^, ll^ to faciliate the derivation. The first 
is the relation (^lexpA^) = ^{exp{Xs + aidi))ai=o- 
In combination with the cumulant expansion theorem 
(expX) = exp{J2{X"')c/n\), we have 

(5iexpA,) = (cxpA,)^I^i^r!^l)£ 

n! 

ri>l 
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(expAa) ^ 



(Al) 



n>0 



The second relation is ((5i(52 exp A5) = gafga^ (exp(AA- + 
ai(5i + a2(52))ai=o,a2=o- Again with the cumulant expan- 
sion theorem, we have 



(Ji(52exp A5) = (expAa) X (A2) 

n>0 



7,1 



j>l ■' n>l 

Putting ah pieces together, we obtain 
((1 + 5i)(l + ^2) cxp A,-) = Df''^(k,){\ + e5(r, K)) 
(A,"((5i+(52))c , k^iW 



1 



E 



7,1 



E 

n>0 



7 1 



(A3) 



oi T7I 



Through the cumulant expansion theorem, we have 



1 + e5(r, fc^) 



exp 



= exp 



E 

a+/3>2 



(A 



(A?)c (-A^) 
a! /3! 



(a + 13)1 
+ 0(«*) 



Terms odd in the power of v vanish due to symmetry 
reason. For example, due to symmetry of w o —v, we 
have {ivi-v2f} = 0, {vf) = (vf). So (wiwi)-(w?v2) = 0. 
Terms of 0(i;'^) do not show in Eq. IA41 Expression 
including terms of order can be obtained following 
Eq. M 

Collecting tersm of the same orders together, we have 



((l + ^i)(l + (52)expA5) 
k 

{6162} + i-^(S2Vl ~ 61V2) + -^{V1V2) 

+CArG(r,fc,) + CG(r,fc,)] . 



Two correction terms Cg and Cmg show up. Cg includes 
all high order corrections, if the density and velocity fields 
are Gaussian. It solely arises from the nonlinear real 
space-redshift space mapping. So we denote it with a 
subscript "G" . It has an exact analytical expression. 



k 



{51S2) 
k^ 



l-^{52Vl - 5iV2) 



(A5) 



i?2 



{5iV2){52Vi) exp(fc^(wiW2)/i7^ 



We notice that this analytical result for Gaussian field 
has been derived by 4^] and shown as their Eq. 32. 

As shown in i jllll Cg can be robustly calculated com- 
bining observations without knowing the underlying cos- 
mology and without introducing extra unknown parame- 
ters. The leading order term is 4-th power in the density 
field, with 



^^^Ul52)+^^-^{{v,-V2){5,+52)) 



p 2 fc4 

+ J^{SlV2){S2Vl) + -j^{viV2f 



(A6) 



To the opposite, C'ng solely arises from non- 
Gaussianities. For this reason, we denote it with the 
subscript "NG" . It is the sum of an infinite series of j-th 
order correlations with j > 3. 



(A7) 



Cng = C'jVGj(r, fcz) 

.k,{5i62{vi-V2))c k'f_{{vi-V2)'^{Si + S2))c 



H 



2iJ2 



In the above equation, we only show the explicit expres- 
sion of CNG,j=3- It is valid even if higher order velocity 
correlation terms are included in (Eq. IA4[) . 

So far the results are exact. To proceed, we have 
to make reasonable approximation on Cng- We argue 
that, due to extra suppression ^ 1 in the deeply 
nonlinear region, higher order non-Gaussian terms are 
expected to be smaller than CNG.j=3- So keeping only 
Cj\!G,j=3 should be reasonably accurate. Nevertheless, 
we will check this approximation in the future, through 
N-body simulations. If needed, CNG,j>4: may be included 
in the calculation. 

In the appendix [B] we will prove that CNG,j=3 is 
equivalent to the term A in 53]. Strictly speaking, its 
Fourier transform Cis!G,j=3{k,u) = A{k,u) in the limit 
that vs 0. 



Appendix B: The Cng,z-A relation 

The term A is derived by [s^ as an additive correction 
to the Kaiser formula. [5^ does not distinguish between 
V5 and V5, so the velocity showing up in the expression 
of A is ve = -I- vg. We consider the hmit vg = 0. 
Re-expressed in our notations, A in (53| is 

Mk,u) = jiiv,-V2)[6,-^)i^52~^)) 

exp(-ik • v)Sy (B1) 
= ikz I {{vi - V2)5i62) exp(-ik • (r)(i^r 



k f 

+ ij^ / ((w2(52ViWi — Wif^i V2W2 + <5if2V2'y2 
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-52Vi\/iVi)} exp(ik ■ (xi - X2))d^Xid^X2-^ 



1 



exp(ik • (xi - X2))d'^xid''x2— . 

Here V is the total volume. Since W2^2ViWi exp(ik-Xi) — 
'Vi[v2S2Vi exp(ik • Xi)] — ikzV2S2Vi exp(ik ■ Xi) and since 
Vi(- • •) integrates to be zero, the term U2^2ViUi in the 
above equation can be replaced by the term —ikzV2S2Vi. 
In total we can do the following replacements in Eq. IBll 

W2(52ViUi -ikzV2S2Vi , 

ViSi\/2V2 — J> ikzViSiV2 , 

51U2V2U2 ikz^Sivl , 
(52WiViUi ~ikz^S2vf , 
i;iViUiV2U2 kl^vjv2 , 



U2V1U1V2U2 kl^vlvi 



(B2) 



Comparing to Eq. IA6[ we prove that, in the limit vg — >■ 
0, CngA^.u) = A{k,u). 

On the other hand, our 4 is not equal to the B term 
in due to differences in the methods and differences 
in approximations made. For example, [53^ sets = in 
their Eq. 18. Inclusion of ^ in our derivation brings 
up new terms such as the term {viV2)'^ in Eq. IA7I In 
future works we will test against N-body simulations to 
compare the two results. 



Appendix C: Modeling r''"{k,u) 



We then have 



P[!{k,u) - Pss{k){l + fW{ky) + C^^'Gik,u){C3) 



The high order correction term 
C]v^G(fc,u) - yxPss{k)fWu^ 
^ (<5iJ2A?). 



E 

n>l 



E 



+ 



Here, y = E,>i(Ai'^i)c/j! = E,>2 (M'^i>c/j!- We can 
resum C]^q in order of the power of S, 



J>3 



(C5) 



{6162X1)0 



(Af52>; 



2! 



The last expression only shows C^q 3 . Finally we obtain 
the expression for r''*, 



1 



Pss(k)(l+fW{k)u'^ 



Pe^B^k)u-i+CNG,3{k,u)+Ca{k.u) 



(C6) 



1 + 



Ps5(k)(l+fW{k)u-^ 



Evaluating r'^'^ requires -83, W and PegBs- The first 
two are observables, as discussed in i jllll The last two can 
be inferred from the observed Pg^ ( §111 Y)\ . So there is no 
uncertainty involved in predicting r''". This will also be 
the case for the window function (Eq. [35]) required 
for velocity reconstruction. 



Following the derivation of Pgg , the real space -redshift 
space density cross power spectrum PJ"| is given by 



PUik.u) = J {{l + 6i)il + 62)eMikzVi./H]) (CI) 
exp(— ik • r)(fir 



(1 + (5i)(l + J2)exp i 



kzVis,, 



H 



exp(-zk-r)d3r^i:)|OG!(^^)£,FOG(fc^) 

Replacing in Eq. IA3I with Ai, we obtain 

((l + ^i)(l + (52)expAi) = (expAi) 
^ ^ y {Xn5i+52))c ^ y {6i62y^)c ,^2) 



{K6i)c {^162)0 



7! ^-^ ?i! 

j>l n>l 



Appendix D: An alternative approach to reconstruct 
the 3D velocity field 

As shown in i jlVl the key to reconstruction the 3D V5 is 
to infer the correct window function . The approach 
discussed in i jlVI is straightforward. But it relies on the 
RSD modeling and is hence susceptible to inaccuracies 
therein. Here we propose an alternative to simultane- 
ously estimate and reconstruct V5. It avoids mod- 
eling the induced RSD, the most difficult part in the 
RSD modeling. So it is less susceptible to uncertainties 
in the RSD modeling. 

The guideline is that, if we construct and V5 cor- 
rectly, we can move particles/galaxies back to their real 
space positions. This will eliminate the V5 induced RSD 
and hence reduce anisotropics in the power spectrum af- 
ter moving (hereafter we denote it as P^^^J^). To fur- 
ther demonstrate this point, let us consider the limit of 
V5 — 0, vb = 0, = and no measurement noise. 
Under this limit, a correct guess of W will faithfully re- 
cover V5. Moving the particles back to their real space 
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positions using this vg, anisotropies in Pinovc(k) will be 
completely eliminated and we recover the isotropic real 
space power spectrum (P„iove(k) = P{k)). This suggests 
that, by tuning until the Pniovc(k) reaches isotropy, 
we can recover the correct and hence recontruct the 
velocity correctly. 

The real situation is more complicated, due to the fact 
that vs 7^ 0, vb 0, Sg ^ and the existence of mea- 
surement noise such as shot noise in the galaxy num- 
ber distribution. However, none of them is correlated 
with the real space density and none of them can cause 
anisotropic pattern the same as . The reconstructs 
V5 without multiplicative error will cause 



P„. 



{k,u) = / (1 + (<5i<52))Anovc(fc.,r)e-^'^-"-rf^r (Dl) 



Like P'*, Pmovo only depends on k and u. So we write 
these dependences explicitly. Here 



We also define a Cmovo through 



1 ~t~ €move(^2 ; ^) 



^move (kz 1 ^) 



,(D2) 



(D3) 



Since is uncorrelated with V5 and v^, we have 

D^rS,?. = Dl^^{kz){e'^^''-s'")\ Here, /^vs = Vz,s - 
g. So far we have neglected shot noise in the galaxy 
distribution. Since it does not correlate with other com- 
ponents, it only causes damping. This effect can be com- 
pletely described by a damping function Dj^f{ku), the 
same as the case of vb ■ For brevity, we will not consider 
this measurement noise hereafter. 

A more complicated issue to deal with is v^, aris- 
ing from the stochastic component in the redshift space 
overdensity (5" (Eq. [55]) . V5 and Vg are correlated, so 
we can not do the average seperately for 

Vz^S and vl g. To see their correlation, let us check 
the limit k ^ 0. Now we have S — ^ zVz/H, the 



starting point to derive the Kaiser formula. Since 5s = 
{5-VzVz,s/H)-VzVz.,s/H,v^eQhi&in5% ~ -VzVz,s/H. 
Through the relation 6'g(k) = (5|(k)l^*(k), we have 
v|(k) ~ (vs(k) • k)u'^W''{'k)k. Notice that the veloc- 
ity field Wg is statistically anisotropic. Under this limit, 
v| is completely correlated with v*. The vs-v| relation 
beyond this limiting case is too complicated to discuss 
here. In future works we will use N-body simulations to 
quantify it. 

Using the relation 

1 + e„ikz,r) = e'liAv,sA..s)/H^+OiA.^s) , (D4) 

we obtain 



.{k,u) = [Pss{k) + PAvAv{k,u)u^ 



xD 



FOG 



(kz) 



(D5) 



Here, PAvAvik, u) is the power spectrum of Avg = V5 — 
v^. Due to the intrinsically anisotropic Vg, PavAvO^) is 
also anisotropic and depends on both k and u. At large 
scale limit, 

PAvAvCk) ~ Pg,e,ik){l-U^W%k))^ 

^ Pososik) ,^^}- ■ (D6) 

The situation beyond this limit is too complicated to dis- 
cuss in analytical ways here and will be postponed for 
future study. Nevertheless, PavAv(}^) has a generic prop- 
erty PAvAv{k,u)/Pss{k) when fc ^> 0, since both V5 
and Vg vanish at large scales. 

We then propose that, by tuning W such that Pmovc 
follows a form like Eq. ID5| we could obtain the correct 
and hence the correct field. So far we still lack of 
a rigorous mathematical proof nor numerical verification 
for the above proposal. However, given its potential in 
reconstructing the 3D peculiar velocity in a less model 
dependent way, we hope to explore this possibility in fu- 
ture works. 
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